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We study, both theoretically and experimentally, modulational instability in optical fibers that 
have a longitudinal evolution of their dispersion in the form of a Dirac delta comb. By means of Flo- 
quet theory, we obtain an exact expression for the position of the gain bands, and we provide simple 
analytical estimates of the gain and of the band widths of those sidebands. An experimental valida¬ 
tion of those results has been realized in several microstructured fibers specifically manufactured for 
that purpose. The dispersion landscape of those fibers is a comb of Gaussian pulses having widths 
much shorter than the period, which therefore approximate the ideal Dirac comb. Experimental 
spontaneous MI spectra recorded under quasi continuous wave excitation are in good agreement 
with the theory and with numerical simulations based on the generalized nonlinear Schrodinger 
equation. 


I. INTRODUCTION 

Modulational instability (MI) refers to a process where 
a weak periodic perturbation of an intense continuous 
wave (GW) grows exponentially as a result of the inter¬ 
play between dispersion and nonlinearity. MI constitutes 
one of the most basic and widespread nonlinear phenom¬ 
ena in physics, and it has been studied extensively in 
several different physical systems like water waves, plas¬ 
mas, and optical devices [i|. For a cubic nonlinearity, 
as the one occurring in the nonlinear Schrodinger equa¬ 
tion used to model optical fibers, the underlying physical 
mechanism can be understood in terms of four-wave mix¬ 
ing between the pump, signal and idler waves. However, 
the scalar four-wave interactions in a homogeneous fiber 
can be phase matched, and hence efficient, only in the 
anomalous group-velocity dispersion (GVD) regime. In 
the normal GVD regme, on the other hand, MI can occur 
in detuned cavities [i, Q , thanks to constructive interfer¬ 
ence between the external driving and the recirculating 
pulse. Alternatively MI with normal GVD can also arise 
in systems with built-in periodic dispersion , among 
which dispersion oscillating fibers (DOFs) have recently 
attracted renewed attention [ilil. In this case, phase 
matching relies on the additional momentum carried by 
the periodic dispersion grating (quasi-phase-matching). 
The occurrence of unstable frequency bands can then be 
explained using the theory of parametric resonance, a 
well-known instability phenomenon which occurs in lin¬ 
earized systems for which at least one parameter is varied 
periodically during the evolution [T^. Up to now, most 
experimental investigations realised in o ptic al fibers have 
been performed with basic sinusoidal |§-[l^ or amplitude 
modulated [l^ modulation formats. In this work, on the 
other hand, we study a radically different periodic mod¬ 
ulation of the GVD, in the form of a periodic train (or 
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comb) of Dirac delta spikes. This is a fundamental and 
widespread modulation format, encountered in a variety 
of physical systems. In optics, delta combs have been ex¬ 
ploited to model lumped amplification in long haul fiber 
optic transmission systems [IS,[HI, oi* to model the power 
extraction in soliton based fiber lasers [T9|. Moreover, 
comb-like dispersion-profiled fibers have been exploited 
to generate trains of solitons starting from a beat sig¬ 
nal [ 2 ^. At more fundamental level kicked systems are 
widely investigated as a paradigm for the emergence of 
chaos in perturbed Hamiltonian systems, with the delta- 
kicked rotor being the most renowned example [T5|. Its 
quantum version is described by a Schrodinger equation 
forced by a Dirac comb and has been extensively analyzed 
to study chaos in quantum systems [T6|. Recirculating 
fiber loops have been used to reproduce the quantum 
kicked rotor with an optical system, to study chaos and 
Anderson localization [^, , and to illustrate how an 

optical system can be used to mimic other physical sys¬ 
tems that are more difficult to reproduce experimentally. 
In the same vein, we hope the experimental setup we 
propose in this paper could be used as an experimental 
platform to investigate such phenomena in the presence 
of nonlinearities, a topic of much current interest. Fi¬ 
nally the approach that we propose to analyse MI in the 
fiber with delta-kicked GVD allows us, on one hand to 
enlighten the featuress of the parametric resonance that 
are not dependent on the specific format of the modu¬ 
lation, and on the other hand to compare and contrast 
the features of the ideal delta-kicked profile with other 
formats including non-ideal (physically realizable) kick¬ 
ing as well as widely employed profiles such as oscillating 
GVD. 

The paper is organized as follows. In Section [II] we 
provide a simple argument allowing to determine the 
central frequencies of the unstable sidebands for general 
periodically modulated fibers. In Section lllli we then 
use Floquet theory to analytically compute the width of 
the gain bands and as well as their maximum gain for 
dispersion-kicked fibers. In Section HVj we investigate nu- 
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merically the effect of the smoothing of the delta comb. 
In Section [Vl we describe the experimental set-up and 
we compare the experimental results with theory and 
numerical simulations based on the generalized nonlin¬ 
ear Schrodinger equation. We draw our conclusions in 
Section ED 


II. IDENTIFYING THE GAIN BAND 
CENTRAL FREQUENCIES 


Consider the NLSE 

.du P 2 {z) d‘^u 


'dz 


■ j{z)\ufu = 0. 


( 1 ) 


We will assume the dispersion I32{z) and the nonlinearity 
coefficient j{z) are of the form 


in the present situation is the two by two real matrix 
defined by (g(0),p(0)) = {q{Z),p{Z)). As a 

result {q{nZ),p{nZ)) = ,^^"(g(0),p(0)). Note that 

iiGcessarily has determinant one, since it is 
obtained by integrating a Hamiltonian dynamics, of 
which we know that it preserves phase space volume. 
As a consequence, if A is one of its eigenvalues, then 
so are both its complex conjugate A* and its inverse 
A“^. This constrains the two eigenvalues of 
considerably: they are either both real, or lie both on 
the unit circle. Now, the dynamics is unstable only if 
there is one eigenvalue A satisfying |A| > 1, in which case 
both eigenvalues are real. We will write A± for the two 
eigenvalues of We are interested in studying 

the gain, that is 

= ln(max{|A+|, |A_|})/Z (5) 


I32{z) = ^av + Pmfziz), ^{z) = 7av + 7m9z{z), (2) 


where fz and gz are periodic functions of period Z 
such that min fz = —1 = min^^, and f^^^^i^fz{z)dz = 

I-z% 9 z{z)dz = 0. 

Let uo{z) = \/Pexp(iP j(z')dz') be a stationary 
solution of m- We consider a perturbation of uo{z) in 
the form u{z^t) = {v{z,t) -h l)uo{z)^ where the pertur¬ 
bation v{z^t) satisfies |u| ^ 1. Inserting this expression 
in o, and retaining only the linear terms we find 


.dv 


% 


dz 


P2{Z) d^V 
2 dt^ 


+ 'y{z)P{v + V*) = 0. 


( 3 ) 


Writing v = q-\-ip^ with q and p real functions, we obtain 
the following linear system: 


dz 
dp 
< dz 


(32{z) d'^p 
2 dt‘^ 


j32{z) d'^q 
2 dP 


— 2^{z)Pq = 0. 


as a function of cj, prn and 7 ^. It measures the growth of 
{q{nZ),p{nZ)). The gain vanishes if the two eigenvalues 
lie on the unit circle. A contour plot of the gain in the 
{uj^Prn) plane, for the case of the delta comb dispersion 
modulation that is the main subject of this paper, can 
be found in Fig. [2l The regions where the gain does 
not vanish are commonly referred to as Arnold tongues. 
We will explain below that, whereas their precise form 
depends on the choice of /z, the position of their tips 
does not. 

Since the system is not autonomous, it cannot be 
solved analytically in general. Nevertheless, the above 
observations will allow us to obtain some information 
about its (in)stability for small 7 ^, and valid for 
all perturbations whatever their specific form. 

To see this, we first consider the case /3m = 0 = 7 ^. It 
is then straightforward to integrate the system (|4]) . The 
linearized Floquet map is then given by 


^2 


UJ 


^lin _ 
^Z,0 — 


(^2 


cos{kZ) —^— sm 

— sin{kZ) cos{kZ) 


{kZ) 


■■= L, ( 6 ) 




Finally, taking the Fourier transform of this system in 
the time variable t, leads to 


dq 

dp 

dz 


2 

I32{z) 


iO^P = 0 , 

uP'q - 2j{z)Pq = 0, 


( 4 ) 


where we used the definiton q{z^uj) = 
J q{z,t)e~^'^^ dt. Note that this is a Hamilto¬ 
nian dynamical system in a two-dimensional phase 
plane with canonical coordinates (g,p). Analyzing 
the linear (in)stability of the stationary solution uo{z) 
therefore reduces to studying the solutions to (j4]) for 
each UJ. Since the coefficients in the equation are 
z-periodic with period Z, Floquet theory applies. This 
amounts to studying the linearized evolution over one 
period Z, to obtain the Floquet map which 


where 

(7) 

Here (32 = /dav > 0 (normal average dispersion), since we 
restrict our investigations to the defocusing NFS. Note 
that the matrix L has determinant equal to 1, as ex¬ 
pected. The eigenvalues of L can be readily computed 
as 


\±{uj,Pm = 0 = 7 ^) = exp(±i/cZ). (8) 

What will happen if we now switch on the interaction 
terms fz{z) and gz{z)l It is then no longer possible, 
in general, to give a simple closed form expression of the 
solution to 0, which is no longer autonomous, and hence 
of the linearized Floquet map Nevertheless, 
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we do know that, for small the eigenvalues of 

7 ^ must be close to the eigenvalues \±{oJ^Prn = 
0 = "fm)- We then have two cases to consider. 

Case 1. k G Z. Now A_(co’,/3^ = 0 = 7 ^) = 

X'^{uj,l3rn = 0 = 7m) 5 they are distinek and they both 
lie on the unit circle, away from the real axis. They 
then must remain on the unit circle under perturbation 
since, for the reasons explained above, they cannot 
move into the complex plane away from the unit circle. 
Consequently, in this case, the stationary solution uo{z) 
is linearly stable under a sufficiently small perturbation 
by fdmfzi^) and and this statement does not 

depend on the precise form of fz{z) or of gz{z). In 
fact, with growing and/or 7 ^, the two eigenvalues 
will move along the unit circle until they meet either at 
— 1 or at +1 for some critical value of the perturbation 
parameters. Only for values of the latter above that 
critical value can the system become unstable. A 
pictorial description of this situation is shown in the left 
hand side of Fig. [H 

Case 2. k = £ e Z. Now A+ = A_ = ±1 is 

a doubly degenerate eigenvalue of Uader a small 

perturbation, the degeneracy can be lifted and two real 
eigenvalues can be created, one greater than one, one less 
than one in absolute value. The system has then become 
unstable! Of course, it will now depend on the type of 
perturbation whether the system becomes unstable, re¬ 
mains marginally stable (the two eigenvalues don’t move 
at all, but stay at 1 or —1), or becomes stable (the two 
eigenvalues move in opposite directions along the unit 
circle). A pictorial description of this situation is shown 
in the right hand side of Fig. [TJ For the Dirac comb mod¬ 
ulation of which is our main object of study in this 

paper, the details are given in the next section. 

In conclusion, examining 0 , one sees that only if 00 = 
uj£, where 

“ £ (\/ (^) ■ 

can an infinitely small Hamiltonian perturbation of 
lead to an unstable linearized dynamics near the fixed 
points uo{z) considered. These values of uj therefore cor¬ 
respond to the tips of the Arnold tongues, that is, to the 
positions of the (centers of) the unstable sidebands of 
the defocusing NFS under a general periodic perturba¬ 
tion fz^gz- This is illustrated for a Dirac comb modu¬ 
lation of the GVD in Fig. O One also observes in that 
figure that, for a value of uj close to some the system 
becomes unstable only for a small but nonzero critical 
value of that we shall compute below for the Dirac 
delta comb GVD. 

Equation ([9|) was derived in [TsI by appealing to the 
theory of parametric resonance and Poincare-Lindstedt 
perturbation theory. Our argument above is elementary 
and shows in a simple manner that the resonant frequen¬ 
cies uj£ do not at all depend on the form of fz or gz- 
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FIG. 1. Sketch illustrating, in the complex plane, the effect 
of the interaction terms fz(z) and gz{z) on the eigenvalues of 
the linearized Floquet map (|6|). Black dots correspond to the 
unperturbed eigenvalues lying on the unit circle (dashed line). 
Coloured dots show the new position of the eigenvalues after 
switching on the perturbations, leading to a stable regime 
when k ^ ^ and an unstable one when k — 


Note that, if fz(z) = sin(^ 2 ;) and gz{z) = 0, a case 
considered in [9|-[i3, the system 0 is equivalent to the 
equation of a harmonic oscillator of (spatial) frequency /c, 
sinusoidally modulated with period Z. In that case the 
system leads to a Mathieu equation for which it is known 
that resonance occurs when the period of the modulation 
is a integer multiple of the half (spatial) period of the os¬ 
cillator, which is 27T/k. 

Additional physical insight can be obtained by ex¬ 
panding Equation for small power, i.e. assuming 
lavP ^ \£\'k/Z. At zero order we recover the well known 
quasi-phase-matching relation Bilil 

+ ‘^lavP — ( 10 ) 

Equation (pT]) entails the conservation of the momentum, 
made possible thanks to the virtual momentum carried 
by the dispersion grating, of the four wave mixing inter¬ 
action between two photons from the pump, going into 
two photons in the symmetric unstable bands at lower 
(Stokes) and higher (antiStokes) frequencies with respect 
to the pump. 


III. CALCULATION OF THE MODULATIONAL 
INSTABILITY GAIN BANDS: DIRAC COMB 


We now turn our attention to the computation of the 
gain G(co’), in particular for values of uj close to the res¬ 
onant frequencies. We concentrate on the special case 
where the GVD is a Dirac delta comb: 


fziz) 


'^5{zlZ-n) 


-1, 


7m = 0. (11) 


Since in the rest of this paper, 7 ^ = 0 , we will drop 
it from the notation. To compute the gain, we need to 
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FIG. 2. (a) Level plot of the gain G{uj,pm) — 

In (max{|A+|, |A_ |}) /Z in the {uj,l3m) plane, for = 1, Z = 
1, 7av = 1 and P = 1. The dashed black lines corresponds 
to the tips of the Arnold tongues Q at a; = cji = 2.1433 
and UJ 2 = 3.2748. The solid red lines corresponds to the gain 
bandwidth, which can be computed from (ITOl) . (b) MI gain 
for f3m = 0.1; Red circles, estimates of maximum gain ([TSl) : 
Black crosses, estimates of the bandwidth (nni. (c) Solid blue 
curve, MI gain for cu = ui; Dashed red curve, approximation 
of maximum gain (1181) . 


[/^av ^m]’) and 

^ I cos - sin 

\sm{pm^Z^ COs[j3rn^z'^ 

The characteristic polynomial of LK is given by 

A^ — 25{uj^ + 1 = 0, 

so that the eigenvalues of m can be computed explicitly 
as: 



A±(w,/3to) = 5{LO,/3m) ± “ 1, (14) 


with 

_L ^ p 

S = cos(kZ) cos (0) - - - - ——— sinfkZ) sin (0 ), 

k 

and 0 = Pm^Z. 

A Taylor expansion of 6{uj,Prn) about (cj^, 0) yields 
6{l3m, 11) ^ (-1)^ [1 + - De{co - wef] , (15) 

where 

D( = &ZV+) (7™-P)" + lV (17) 

The dependence in (not in Pm) entails that the sign 
of the kick has no incidence in this regime, i.e. assuming 

Formula ([T5]) shows that (0,cc;^) is a saddle point for 
5{pm^^)- If ^ is even, A+(y5^,cc;) > 1 occurs close to 
( 0 , 0 ;^), and if ^ is odd, X-{Pm,oj) < —1 close to (0,^;^). 
More precisely, 

max(|A+|, |A_|) = 1 + 7^/3^ - De{co - uje)"^ 

from which we can find an estimate of the gain amplitude 
GiuJi^Pm) and of the bandwidth B{uj£,pm) near the tips 
of the tongue at as: 

(18) 


compute the linearized dynamics and determine the 
behaviour of its eigenvalues X±{Pm^^) in the neighbour¬ 
hood of Pm = 0 and uj = uj£ in the (cj, /3^)-plane. 

In this case the linearized Floquet map is easily seen 
to be explicitly given by 

= LK (12) 

where L is defined by Equation (|6]), but now with P 2 = 


B(uj£ , Pm) 


\Pm\ 'JayP 

^ Pf) +havPy 


(19) 


Note that the threshold value for Pm above which insta¬ 
bility occurs can be read off from the above by setting 
\5{pm^^)\ = 1 which corresponds to 


l/3m| > 



Ui\. 
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FIG. 3. MI gain for a; = as a function of Red circles, 
estimated gain given by (ES- Parameter values are /3av — 

Z = 1, 7av = 1, R = 1, Pm = 0.1. 


This confirms again, as expected, that an arbitrary small 
Pm will generate instability right at cj = In Fig. 
[2^ we show an example of the analytically computed MI 
gain, showing the first two Arnold tongues. As can be 
seen, for a small enough strength of perturbation, let’s 
say \pm\ < 0.1, the approximation ([19]) gives a good es¬ 
timate of the width of the parametric resonance (see red 
curves). This situation is detailed further in Figs. [2]3,c, 
showing a section for Pm = 0.1 and cj = cji, respectively. 

Finally, a straightforward calculation gives the asymp¬ 
totic behaviour of the gain G at for I large and Pm 
fixed, that is 



FIG. 4. Gain G(ijJi) G((jJi,Pm), with Pm — 0.1, calculated 
numerically at the first parametric resonant frequency uji , for 
two approximations of delta functions as a function of their 
inverse height h~^. Solid blue curve, gaussian approximation; 
Dashed red curve, rectangular pulse approximation; Dashed 
horizontal line, Dirac delta limit; Dash-dotted horizontal line, 
sinusoidal modulation. Parameter values are /3av = 1, Z = 1, 
7=1, P=l. 


minimal value, and even the fact that an instability is 
indeed generated, does depend on the precise form of fz- 
For the Dirac comb the explicit expression for the gain 
in this regime can be read off from ([2Q|) . 

IV. APPROXIMATIONS OF THE DELTA 
FUNCTION 


^ , J 4/3^sin^(a(^))(7avP)2 - 

^ - 

( 20 ) 

with a{i) = whenever AP‘^^ sin^(a(^)) — 

> 0, and G(uJi) ^ 0 otherwise. 

In Fig. [3l we show an example of the analytically com¬ 
puted MI gain at uji as a function of 1. We compare it 
to the approximation ([20]) . which is very accurate, even 
for small £ (see red circles). Note in particular that the 
oscillating behaviour of the gain is well captured by (|2Q|) 
which, for £ large enough and pm small, can be approxi¬ 
mated by 


G{uJi,Pm) — \Pn 




^tt£ 

Pav 


In order to shed light on the dependence of the gain 
on the shape of the periodic modulation, and also with 
an eye towards the experimental realization of the Dirac 
comb fiber, we now analyze what happens when the Dirac 
comb is approximated by a train of physically realizable 
“kicks”. We thus consider a smoothened Dirac comb de¬ 
scribed by 


fz{z) 



z — nZ) 


- 1 , 


( 21 ) 


where we normalize the positive function rj{z) in order 
to have fz{z)dz = 0. For a rectangular pulse of 

width re < Z, we get 


7]{z) 



-f <2< f 
elsewhere. 


( 22 ) 


Summing up. It is clear from the above that, precisely 
at the values which only depend on Z and on ^avP^ 
but not on the precise form of fz , any small perturbation 
can create an instability and hence a gain. At frequen¬ 
cies uj near these particular values, a minimal threshold 
strength of Pm is needed to create an instability. This 


For a gaussian function r]{z) = hexp(—2)^/(2rc^)), the 
maximum amplitude of the kick can be calculated as 
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that in the limit w Z gives 



Note that, in these models, we have {/3m > 0) 


Pav — /^min ^ /^av 5 /^max — /^av 3~ {h l)/5?7 


Hence h = 2 corresponds to a rather symmetric situ¬ 
ation where ydav is close to the midpoint between /3min 
and ydmax, so that p 2 {^) fluctuates symmetrically about 
its average value. Whereas h >> 1 corresponds to a 
very asymmetric situation where P 2 {z) has a large abrupt 
peak. The parameter h therefore controls the shape of 
the GVD modulation at fixed and /3min (or Pm)- 
As shown in the previous section, by changing the 
shape of the kick, we do not change the frequency of 
the parametric resonances. The smoothing of the delta 
function nevertheless does modify the characteristics of 
the MI by changing the value of the gain, as we now 
illustrate by computing the gain numerically at the res¬ 
onant frequencies uj£. An example of how the changing 
shape of the modulation fz modifies the first parametric 
resonance is illustrated in Fig. 01 that shows the gain 
G{uJi,Pm) at fixed uji and as a function of the peak 
amplitude h (or, equivalently, the width w) of the kicks. 
We make the following observations. First, a good ap¬ 
proximation of the gain given by the Dirac comb is ob¬ 
tained for h > 10, both for the rectangular and gaussian 
pulses. Second, for h = 1, the gain of the square pulse 
modulation is zero, as expected, since we are then in the 
limit case of a constant modulation (and normal GVD). 
Third, it is apparent that the Dirac comb gives the high¬ 
est possible gain, for a fixed area of the kicks and fixed 
Pm and ydav Finally, it is interesting to note that a sinu¬ 
soidal modulation /z(^) = sin(^ 2 ;), with the same value 
of Pm^ gives a gain close to one half with respect to the 
delta case. Indeed for a sinusoidal modulation, it has 
been shown that (see Equation (7) from Ref. [l^ ) 


Gsin{^£’i Pm) — 


A 


' 2 iri 


lavP- 


(23) 


By expanding Equation (|23|) for small Pm^ we get 


Gs\n{^£’)Pm) — 2^G{uj£^Pm)^ 


(24) 


at first order in Pm- In conclusion, a large concentrated 
perturbation of the GVD about its average enhances the 
MI gain. 


V. EXPERIMENTAL RESULTS 

It is well known that in homogeneous fibers, the GVD 
depends on the diameter of the fiber. One can therefore 
modulate the GVD by modulating the diameter of the 
fibers as a function of z, as in We manufactured 


three different microstructured optical fibers modulated 
by a series of Gaussian pulses to approximate the ideal 
Dirac delta comb studied in Section imi The change of 
their outer diameters d{z) along the fiber is represented 
in Eig. [5fa). As can be seen in the inset, their diameters 
have a gaussian shape with a standard deviation re, which 
is the same for all three fibers, and very small {w 0.14 
m) compared to the period of the comb (10 m). Hence 
we can write 


d{z) — dmin T 3g^{z) ^ dmax — dmin “t“ ^5 

where gn,{z) = y„exp(-i The three fibers 

have a very similar minimum diameter dmin = 137/im, 
while their maximum values are different. We have 
^max —172 /im, d^ax —207 jj^m and d^ax —240 gm for 
fibers labelled A, B and G, respectively, corresponding 
to 5^ = 70/rm, 5^ = 103/im. To understand 

how the two experimental parameters w and S control 
the quality of the approximation of the delta function on 
the one hand, and the value of Pm on the other hand, 
we proceed as follows. Eirst, dav = do + so 

that d^ = 139.1, d^ = 141.2, d^ = 143.2. A first order 
Taylor expansion of /32(d) about dav yields 


182(2) = j82(d(z)) ~ ,82(da,y)+,82(da,y)S g^u(z) - 


Gomparing this to (|2j) we find 


6w 


Pav — P2{dav)<, Pm — P 2 {^aiv) ry V^2^, 


and 


fz(z) = — J=gw(z) - 1. 
rev 27r 


Hence, with the notation of Section llVl h = —28.5. 

This corresponds to ^ 0.03 proving that these Gaus¬ 
sian pulses should induce a very similar parametric gain 
compared to ideal Dirac delta functions (see Eig. El). 
Eurthermore, the height 6 of the Gaussian pulse controls 
Pm, which will allow us to investigate the impact of this 
parameter on the first MI side lobe gain, as it was done 
in the theoretical study and illustrated in Eig 121 Eor all 
three fibers, the ratio of the diameter of the holes over A 
(the pitch of the periodic cladding) is assumed to be con¬ 
stant along the fiber and estimated to about 0.48 from 
scanning electron microscope images. The diameter vari¬ 
ations of the fibers are proportional to those of the pitch, 
with Amin = 3.445 /rm corresponding to the minimum 
value of the diameter (137/im, green line in Eig. [5fa)) 
and AAP^ 4.48, 5.17 and 6.03 (im, blue, black and 
red curves respectively, for the maximum values. As an 
example, the Group Velocity Dispersion (GVD) curve 
corresponding to the minimum pitch value has been cal¬ 
culated from Ref. [23 and is represented in Eig. [Sfb) as 
a green curve. Its Zero Dispersion Wavelength (ZDW) is 
located at 1055 nm while those of the GVD curves cor¬ 
responding to the maximum values of the diameters of 
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FIG. 5. (a) Outer diameter evolution versus fiber length for 
fibers A, B and C. Inset: Snapshot on one leitmotiv, (b) 
Calculated GVD curves corresponding to the background di¬ 
ameter (green curve) and the maximum diameter values of 
fibers A (blue curve), B (black curve) and C (red curve), (c) 
GVD and nonlinear coefficients evolution versus the diameter 
at 1052.5 nm. Curves in (b), (c) and in (d) calculated from 
Ref. [ 2 ^ with d/A = 0.48 and Amzn = 3.445/im. 


fibers A, B and C are red-shifted to 1110 nm, 1136 nm 
and 1168 nm respectively (Fig. Efb)). In order to give 
another illustration of the large dispersion variations in¬ 
duced in these fibers by varying their diameters, the max¬ 
imum GVD values for fibers A, B and C have been calcu¬ 
lated at a fixed wavelength (1052.5 nm) and compared to 
the background value. As can be seen in Fig. [5fb), an in¬ 
crease of the diameters by a factor of only 1.27, 1.45 and 
1.75, leads to a one order of magnitude improvement on 
the GVD values: 21, 29 and 35 respectively. Under such 
large variation of the fiber diameter, the GVD can no 
longer be considered as proportional to the pitch value, 



FIG. 6. Scheme of the experimental setup. MOD: intensity 
modulator; YD FA: ytterbium doped fiber amplifier; DOF: 
dispersion oscillating fiber; PM: polarization maintaining; 
TL: tunable laser; OS A: optical spectrum analyzer. 


as was the case in Ref.[9[ for instance. This is illustrated 
in FigISfc), where the evolution of the GVD calculated 
at 1052.5 nm as a function of the fiber diameter is rep¬ 
resented. It can be seen to be well approximated by 
an affine function in the range between 140//m and 180 
/im, but not beyond. As a consequence, the shape of 
the GVD variations will be slightly different from the 
one of the diameter, specifically for fiber C. However, we 
checked numerically that this can be considered as rela¬ 
tively weak distortions that do not significantly impact 
the gain of the MI process. We can still consider that the 
key parameters remain the different heights in GVD of 
the Gaussian-like pulses in fibers A, B and C. In order to 
get a more complete picture of the impact of the fiber di¬ 
ameter variations on its guiding properties, the variation 
of the nonlinear coefficient is plotted as a function of the 
fiber diameter in Fig. [5fd). The most important feature 
to note here is that the amplitude of variation is much 
smaller, and only the same order as the one of the diame¬ 
ter itself. Hence, these variations are more than one order 
of magnitude lower than those of the GVD and we have 
checked numerically that their impact on the MI process 
is negligible. Consequently, we can infer that these fibers 
represent a good prototype to validate our theoretical in¬ 
vestigation in the previous sections, where only longitu¬ 
dinal GVD variations have been taken into account. The 
experimental setup is schematized in Fig. [6l The pump 
system is made of a continuous-wave tunable laser (TL) 
diode that is sent into an intensity modulator (MOD) in 
order to shape 2 ns square pulses at 1 MHz repetition 
rate. They are amplified by two Ytterbium-doped fiber 
amplifiers (YDFAs) at the output of which two succes¬ 
sive tunable filters are inserted to remove the amplified 
spontaneous emission in excess around the pump. These 
quasi-CW laser pulses have been launched along the bire- 
fringent axis of the fibers. The pump peak power has 
been fixed to 6.5 W and the pump wavelength at 1052.5 
nm for fiber A. The output spectrum recorded at its out¬ 
put is represented as a blue curve in Fig. [3 (a). Two MI 
side lobes, located at +/-4.8 THz appear on both sides 
of the pump. These experimental results have been com¬ 
pared with numerical simulations performed by integrat¬ 
ing the generalized nonlinear Schrodinger equation. We 
used the GVD, ^^4 and 7 variations calculated from the 
measured diameter values (see Figs. 0 . Other parame¬ 
ters are extracted from experiments and are listed in the 
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caption of Fig. [3 Note that we checked that except the 
longitudinal GVD variations, all other parameters can 
be assumed to be constant and equal to the average val¬ 
ues. As can be seen in the blue curve in Fig. El^b), two 
symmetric MI side lobes also appear in the simulated 
spectrum, in a very good agreement with experiments. 
Their positions have been compared with the predictions 
of Equation [HI represented by green dashed lines in Fig. 
[3(b) (calculated with 59 ps‘^/km). An ex¬ 

cellent agreement is also obtained. In order to show that 
the MI gain is larger when the weight of the Dirac delta 
function is increased, we performed similar experiments 
in fibers B and C where the areas of the Gaussian pulses 
are larger than in fiber A. However, in experiments, due 
to the fact that the Dirac comb has been approximated 
with a series of Gaussian functions, changing their am¬ 
plitudes also modifies the average value of the dispersion. 
As a consequence, MI side lobes would be generated at 
different frequency shifts. In order to keep constant the 
position of the MI side lobes, and then provide a cor¬ 
rect comparison with the theoretical study, one have to 
take care to keep the average GVD value constant in all 
the fibers. To do so experimentally in fibers B and G, 
we slightly tuned the pump wavelength until the first MI 
side lobe is generated at 4.8 THz as in fiber A. As can 
be seen in Fig. [3(a) (red and black curves), the position 
of the first MI side lobe in fibers B and G is indeed lo¬ 
cated at that frequency by tuning the pump wavelength 
to 1061.8 nm and 1067 nm, respectively. We can there¬ 
fore consider that the average GVD values are very sim¬ 
ilar in the three fibers and hence that only the areas of 
the Gaussian functions, i.e. the equivalent of the Dirac 
weights, vary. The amplitudes of the first MI side lobes 
generated in fibers B and G are indeed larger compared to 
fiber A, as predicted by the theory. This is in pretty good 
agreement with numerical simulations (Figl3(b)), where 
the same procedure was used. We found that the average 
values in fibers B and G are _ _|_q 53 ps‘^/km 

and _ _^Q 52 ps^/km respectively. The small 

discrepancy between these values is attributed to spuri¬ 
ous longitudinal fluctuations arising during the drawing 
process. Indeed, as can be seen in Fig. [SKa), the back¬ 
ground over which Gaussian pulses are superimposed is 
not perfectly fiat, and in fibre G, it is not horizontal. To 
counterbalance these imperfections, it was necessary to 
adjust the average dispersion values. We checked that 
with a series of perfect Gaussian pulses superimposed on 
a fiat and horizontal background, the same average GVD 
value would be obtained. Furthermore, we can note that 
in fibers B and G, additional MI side lobes are generated 
due to the periodic modulation of the GVD (labelled Mli 
in [3(a), up to 5 in fiber B). Their positions are also well 
predicted by numerical simulations and by Equation [9] 
(green lines in [3(b)). This excellent agreement confirms 
that their positions indeed scale approximately as \//, I 
being the side lobe order, that is the typical signature of 
the MI process occurring in dispersion oscillating fibers. 
It was already reported experimentally in Refs. im, 
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FIG. 7. Spectra recorded at the output of fibers A, B and 
C. We used the same color label than in Fig. [51(a) Exper¬ 
iments and (b) numerical simulations. Parameters: = 

5.6 X /34 = —1 x o = bdB/km. 


with a sinusoidal variation of the GVD, modulated in 
amplitude or not and it is now illustrated in this paper 
with a Dirac delta comb. Moreover, in fibers B and G 
two symmetric side lobes that are not predicted by the 
theory appear around the pump at about 2.15 THz (la¬ 
belled spurious side lobes in EigI3(a)). They result from 
a non-phase matched four-wave mixing process involv¬ 
ing the pump, the first and second MI side lobes. The 
energy conservation relation involving these waves pre¬ 
dicts a frequency shift of 2.2 THz from the pump for the 
fourth wave, that is in good agreement with the shift of 
2.15 THz measured experimentally. 


VI. CONCLUSIONS 

Modulation instability has been investigated theoret¬ 
ically and experimentally in dispersion kicked optical 
fibers. An analytical expression of the parametric gain 
has been obtained allowing to predict the behavior of 
the MI process in such fibers. Specifically, it was shown 
that increasing the weights of the Dirac functions leads 
to larger MI gains for the first MI side lobe. We exploit 
the fact that the Dirac delta comb can be well approx¬ 
imated by a series of short Gaussian pulses in order to 
perform an experimental investigation using microstruc- 
tured optical fibers. We then experimentally report, for 
the first time to our knowledge, multiple MI side lobes at 
the output of these dispersion kicked optical fibers. We 
demonstrate that they originate from the periodic varia¬ 
tions of the dispersion. We also validate experimentally 
that increasing the height of the modulation leads to a 
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larger gain for the first MI side lobe. This illustrates that 
optical fibers constitute an interesting platform to real¬ 
ize experimental investigations of fundamental physical 
phenomena. 
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